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Self-consistent solutions of the Ginzburg-Landau system of equations, which describe the behavior of 
the order parameter i/j, and the magnetic field distribution B, in a long superconducting cylinder of 
finite radius R, in external magnetic field H , when a vortex line, carrying m flux quanta, is situated 
on the cylinder axis (a giant m- vortex state), are studied by using numerical method. If the field H 
exceeds some critical value Hs , the giant m- vortex state solution becomes unstable (in response to 
small fluctuations in its shape) and the solution passes to a new stable edge-suppressed form. The 
quantum number m in this state does not change, but the order parameter diminishes by a jump 
(almost to zero) near the cylinder surface; however, superconductivity remains in the deep, at some 
distance from the cylinder axis. This edge-suppressed state exist in the fields Hi™ 1 * 1 < H < H c m \ 
where He™* is the field, in which the second order phase transition into the normal state occures. 
(In the case of large radii R and finite m the field He coincides with the critical field H C 2, in 
which the superconductivity vanishes in the bulk.) If the magnetic quantum number m is large, 
the edge-suppressed state degenerates into the usual state of surface superconductivity, and may 
survive up to the field H c z = 1.69ff C 2. The magnetic moment, the total magnetic flux, the Gibbs free 
energy and other characteristics of the system are found, as functions of the field H and temperature 
T, for different radii R, quantum numbers m and parameters X of the Ginzburg-Landau theory. 
The intervals of R, T, X, m are found, where the edge-suppressed solutions exist. The giant m- 
vortex states exist in both type-I and type-II superconductors, but the edge-suppressed solutions are 
possible only in type-II superconductors (with K > l/\/2). The paramagnetic effect in mesoscopic 
samples and also the possible connection of the theory and experiment are shortly discussed. 



I. INTRODUCTION 

The behavior of finite size superconductors in mag- 
netic field was studied on the base of nonlinear system of 
Ginzburg-Landau equations [1] in numerous theoretical 
papers (see, for instance, [2-15]). The results of [2-15] 
were used, in particular, to explain some anomalies, ob- 
served in a number of experiments with small size super- 
conductors, placed in external magnetic field [16-24]. 

In our recent paper [25], by using numerical meth- 
ods, the self-consistent solutions of nonlinear Ginzburg- 
Landau equations were studied, for a long superconduct- 
ing cylinder of finite radius in external axial magnetic 
field, with no vortices inside the cylinder (m = 0). It 
was found, that type-II superconducting cylinder, which 
is in the Meissner (vortex-free) state, upon increasing the 
field H, may pass by a jump into a special (also vortex- 
free) edge-suppressed state. The order parameter ip in 
this state is strongly suppressed (practically to zero) in 
the vicinity of the cylinder surface, however, the super- 
conductivity survives near the specimen center, where tp 
remains finite. Such edge-suppressed solutions were not 
studied previously [1-15]. 

Apart from the the Meissner-state (m = 0) more 
general one-dimensional solutions exist, when the vor- 
tex line, carrying arbitrary number of flux quanta (m = 



1, 2, 3, . . .), is situated on the cylinder axis. Such solu- 
tions (which are called "the giant-vortex states") have 
been extensively studied in the literature ([2-30]) and 
used in the interpretation of experiments. The shape 
of the giant-vortex changes smoothly with the external 
field H. As is shown in the present paper, for a cylin- 
der of finite radius, the giant-vortex state (with vortic- 
ity m > 0) becomes unstable, when the field H exceeds 
some critical value i?i . After that, the giant-vortex 
solution transforms its shape by a jump and acquires the 
"edge-suppressed" form with the same vorticity m. Such 
edge-suppressed states are characterized by a strongly 
suppressed value of the order parameter ip near the cylin- 
der boundary, and are analogous to the edge-suppressed 
(or, rim-suppressed) states, found in [25] for the case 
m = 0. Previously, such edge-suppressed modifications of 
the giant-vortex states were not studied in detail, though 
some evidence of the edge-suppression was noticed also 
in [2,3,7]. 

We should stress, that the edge-suppressed state is 
noting else, but the usual giant m-vortex solution in a 
field H > Hs ■ If the field H is increased further, the 
edge-suppressed solution continues to vary its form, the 
maximal value of the order parameter ipl^ x gradually di- 
minishes, and for ip^max *C 1 the edge-suppressed solution 
may be expressed in terms of the Kummer functions (see, 
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for instance, [6,8]). 

As is shown below, the edge-suppressed states can be 
formed only in type-II superconductors, with k > l/\/2, 
where k is the Ginzburg-Landau parameter. In type- 
I superconducting cylinder, with x < 1/V2, the edge- 
suppressed states do not exist. 

It is known, that for m > the order parameter van- 
ishes in the vicinity of the cylinder axis according to the 
law ^(r) ~ r m [29]; if the field increases, -0 vanishes 
also near the surface (where the edge-suppressed state 
forms), so, the superconductivity persists only at some 
distance from the cylinder axis. If the magnetic quan- 
tum number m is increased, the region, where the order 
parameter ip(r) ^ 0, shifts more and more toward the 
cylinder surface. For a fixed cylinder radius (R =const, 
m > 1/V2) there exists some maximal value m — m maxi 
for which the edge-suppressed state degenerates into the 
usual state of the surface superconductivity [26-30] . 

The transition of a superconducting cylinder into the 
edge-suppressed state may be accompanied by jumps in 
the dependencies versus H of such quantities, as: the 
specimen magnetic moment (— 4ttM); the total magnetic 
flux, captured inside the system ($1); the system Gibbs 
free energy (G); and by peculiarities in the behavior of 
other parameters, which characterize the superconduct- 
ing state. These topics are discussed below in detail. 
No immediate comparison of the theory with the exper- 
iment is possible, because the model case of an infinitely 
long cylinder approximates the real samples geometry 
[16-24] very remotely. However, a number of qualita- 
tive predictions, which follow from this model, may be 
used for interpreting some of the peculiarities, observed 
in the experiments with mesoscopic samples. In particu- 
lar, the controversial paramagnetic Meissner effect in su- 
perconductors is shortly discussed below, basing on the 
Ginzburg-Landau theory approach. 

The paper is organized as follows. In Sec. 2, the ba- 
sic equations and the boundary conditions of the problem 
are written. The numeric algorithm, used to find the self- 
consistent solutions of nonlinear system of equations, is 
shortly described. Sec. 3 contains the results of numer- 
ical calculations. We are forced to present our results 
by a large number of graphics, which illustrate different 
details of the system behavior in this many parameter 
problem (the solutions depend on m, x, R, H, the spec- 
imen temperature T, the critical temperature T c and the 
coherence length £o)- The necessary comments accom- 
pany the presentation of the material. Sec. 4 contains a 
short resume and discussion of the results. 



II. THE SETTING OF THE PROBLEM 



cylinder element. The basic system of the Ginzburg- 
Landau equations [1] is of the form 



rotrotA=— j s , — ^ = ^ (^ Ve ~ A 



<Po 



where A is the vector-potential of the magnetic field 
(B = rot A), j s is the current density inside the supercon- 
ductor, A is the field penetration depth, £ is the coherence 
length, A = x£. The order parameter, in a general case, 
is written as ^ = ?/>e ie , where tp is the modulus and 
is the phase of the order parameter. From the singlc- 
valuedness of ^ the condition follows 



V6rfl = 27TTO, 



where the contour C embraces the vortex axis, m is in- 
teger (the topological invariant, or fluxoid, or vorticity), 
which shows, how many vortices may be present inside 
the contour R. 

In the cylidrical system of co-ordinates r,<p,z, with z 
axis directed along the cylinder element [when the vector- 
potential has only one component, A = e v A(r)], these 
equations may be written in the dimensionless form 



d 2 U 
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IdU 
P dp 



-tp 2 U = 0, 
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Here, instead of the dimensioned potential A, field B 
and current j s , the dimensionless quantities U(p), b(p) 
and j(p) are introduced: 
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(The field B in (3) is normalized by H x = <j> / {2n A 2 ) , 
with b = B/Hy. instead of H\ one can normalize by 
= ?W(27r£; 2 ) = H c2 [28-30], or by thermodinamical 
critical field H c = xH\/V%. The coefficents in (l)-(3) 
would change accordingly.) 

The magnetic flux, confined inside the contour of the 
radius r, is 



Consider long superconducting cylinder of radius R in 
the external magnetic field H, which is parallel to the 
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Thus, the potential U(p), in the normalization adopted 
above, is related to the flux &(p) by a simple formula 
(f) = $/0 o = U( P ) + TO. 

Because the magnetic flux through the vanishing con- 
tour is zero, and the field B\ r= R = H, the following 
boundary conditions correspond to Eq. (1): 



U\ 



p=0 



-TO, 



dU_ 

dp 



where Pl = R/X, h x = H/H x , H x = </> /(27rA 2 ). 

As to the Eq. (2), we shall take the usual boundary 
condition on the external surface [1]: d^p/dp\ p=Pl = 0. 
The order parameter at the center is either maximal (if 
to = 0), or zero (if m > 0, see, for instance, [29]), thus, 
the following boundary conditions correspond to Eq. (2): 



dip 
dp 



= 0, 



p=0 



V>lp=o = o, 
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dp 



= (to = 0), 



p=pi 



= (m>0). 
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The magnetic moment (or, magnetization) of the cylin- 
der, related to the unity volume, is 
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where B av is the mean field value inside the supercon- 
ductor, $i = <f>(i?), S — irR 2 . In normalization (3), 
denoting b = B av /H x , h\ = H/Hx, M\ = M/H x , one 
finds from (7): 

- - 2 $i 

4nM x = b-h x , b=^(Ui+m), fa = — = Ui + m, 
Pi 9o 



U 1 = U(p 1 ), Pl = 
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For the difference of Gibbs free energies of the system 
in superconducting and normal states, AG = G s — G n , 
it is convenient to use the exact expression: 

AG = F s0 -±MH+^-m(B(0)-H), (9) 



Hi f 



dv, H c = 
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where corresponds to the superconductor condensa- 
tion energy, M is the cylinder magnetic moment, Z is 
the cylinder length, to is the magnetic quantum number, 
B(0) is the magnetic field at the cylinder axis, H is the 



external magnetic field. The expression (9) follows from 
the general formula, obtained in [31] for the free energy 
of a hollow superconducting cylinder, if the radius of a 
hollow is put to zero (see also [32,33]). Using normaliza- 
tion (3), and also (8), one finds from (9) the normalized 
expression 
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The expressions (3), (8), (10) will be used in Sec. 3 for 
calculating the corresponding quantities. 

We remind, that the field penetration depth A and 
the coherence length £ = A/x depend on temperature. 
Thus, the expressions above depend implicitly on tem- 
perature, and, formally, are valid for arbitrary values of 
T. (Though, the Ginzburg-Landau equations themselves 
are applicable only in the limit T — > T c [34,28-30] , when 
the expression £(T) = £o/\/l — T/T c may be used for 
the coherence length). 

It is appropriate to make here a comment on the itera- 
tion procedure, we used to obtain the self-consistent solu- 
tion of the system of equations (l)-(6). First, some trial 
function ip(p) was chosen, and the solution for U(p) was 
found from Eqs. (1), (4). Then, this function U(p) was 
introduced into Eq. (2), which was solved with account of 
the boundary conditions (5), (6) and new function tp(p) 
was found. Further, Eq.(l) was solved again, and all the 
procedure was repetcd, until the functions ip(p) and U(p) 
ceased to change, representing the self-consistent solution 
of the system (see also [25,35]). However, the numeric al- 
gorithm, used in our calculations, docs not allow to find 
the unstable "superheated" and "supercooled" states in 
type-I superconductors (with k < 1/ \/2) , where other 
methods are required (see [32,36,37]). 

The question of stability requires an additional illuci- 
dation. Alongside with the centrally symmetric giant- 
vortex states, the so-called multi-vortex solutions of the 
same vorticity to exist. The multi- vortex states corre- 
spond to the case, when a single giant vortex (with vor- 
ticity to) decays into a number of vortices (with vortic- 
ities TOi, X) m « = m )- These vortices can be arbitrary 
positioned on the cylinder cross-section, so the multi- 
vortex state may be no longer centrally symmetric. The 
question, which of these multi- vortex states possesses the 
minimal free energy, is important, because the physical 
system is usually assumed to occupy the state with the 
smallest energy (the equilibrium ground state), and all 
other states are considered as energetically unstable (or 
metastablc). The free energy calculations to find the 
equilibrium multi- vortex configuration is a difficult prob- 
lem, which requires the specific methods of investigation. 
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(See, for instance, the theory [30], where the triangu- 
lar lattice of single vortices is found as an equilibrium 
state for the case of infinite superconductor, placed in the 
field H ~ H c i. The numeric calculations for mesoscopic 
suprconducting discs may be found in [6-15], where a 
simplified model of the vortex distribution is considered. 
To find the equilibrium vortex configuration in a general 
case of a finite-size superconductor, placed in arbitrary 
magnetic field H , is a problem as yet unsolved.) 

In the present work a limited task was set: to de- 
scribe the formal properties of the giant-vortex states, 
which were not reported previously (in particular, the 
edge-suppressed transitions and the jumps in the mag- 
netization curves with a fixed vorticity to). The gen- 
eral problem of energetic stability and of the equilib- 
rium vortex configuration is left outside the scope of 
the present investigation. Some of the edge-suppressed 
states, which are described below, might in reality be 
energetically metastable and decay into an equilibrium 
ground state configuration. However, as is well known, 
the physical system not necessarily must be in equi- 
librium, but (with some probability) may also occupy 
the excited metastable states. The giant-vortex and the 
edge-suppressed solutions provide the system this oppor- 
tunity, so the behavior of a superconductor in such states 
might be of some interest. 

The results of our calculations are presented below. 

III. THE NUMERICAL RESULTS 

The results of the investigation are illustrated below 
by a number of graphics, which detalize various aspects 
of this multy-parameter task. (Some of these results can 
be obtained by the methods, used in the papers, cited 
above. However, for the uniformity of presentation, we 
prefer to reproduce them by our method.) 

Evidently, the superconductivity would be destroyed 
by a large external field H with subsequent transition 
to the normal state. In Fig. 1(a) the phase diagram is 
depicted (for x = 2 and different to), which separates 
the normal (n) and superconducting (s) states on the 
plane of variables R\ = R/X and fi£ = H/H^. The phase 
diagram is a many-sheet surface, each sheet corresponds 
to the fixed value of the vorticity to. The sheet to = 
is the lowest laying (according to the increasing to), 
the sheet to = 1 lays above it, still more above lays the 
sheet to = 2, and so on. (It is sufficient to consider 
only the values to > 0, because the results for to < arc 
easily reproducible with account of the evident symmetry 
(m,H) ^=> (-m,-H).) 

On each of TO-sheets of the phase surface there is a 
phase curve (h c , which may be marked by the index 
to), which separates the region of values (R\, h^), where 
the superconducting solutions exist (i.e. V > 0), from 
the region of normal state solutions (tp = 0). For a 



bulk superconductor (R\ 3> 1) the superconducting state 
boundary (in the case of type-II superconductor) lays at 
h ( = H/H<: = 1, i.e. at H = H c2 = 0o/(2< 2 ) [28-30] 
(H C 2 = iff). As can be seen from Fig. 1(a), for small 
R\ the superconducting states (with different values of 
to) are also possible in the fields, exceeding the critical 
value H C 2 (i-c. for > 1). For fixed R\ and to there 
exists a maximal field fi£ = he , above which the super- 
conducting solution is impossible. As is well known, the 
maximal field, at which the superconductivity is still pos- 
sible (for a bulk specimen with R\ » 1) is h^ = 1.69 (i.e. 
H = H c3 , where H c3 = 1.69i? C 2 is the field, at which 
the surface superconduvtivity nucleates in macroscopic 
specimens [26-30]). As is evident from Fig. 1(a), for 
small R\, the superconductivity is also possible in the 
fields, exceeding fi£ — 1.69. [The analogous phase dia- 
gram in the case of type-I superconductor (with x. = 0.5, 
to = 0, 1, 2) is depicted in Fig. 1(&).] 

Inside each of the superconducting regions in Fig. 1(a) 
there exists a singular line h s , which is depicted, for clar- 
ity, on a separate phase m-sheet (see Fig. 2). In Fig. 
2(a) the phase sheet to = (for x = 2) is depicted. The 
solid line h c corresponds to the superconducting state 
boundary. The letter s marks the region, where the su- 
perconducting state exists (ip ^ 0), the letter n denotes 
the normal state (ip = 0). 

To every point (R\, h^), which lays inside the supercon- 
ducting region in Fig. 2(a), corresponds some solution of 
Eqs. (l)-(6). When this point is shifted, the solution 
changes accordingly. The central part of the drawing 
(which lays in the withinity of ~ 0) corresponds to 
the superconducting Meissner state. The order param- 
eter here is ip ~ 1, the external field H is almost com- 
pletely screened out and penetrates only to the distance 
<~ A from the superconductor surface. The shape of the 
Meissner solution depends on H only slightly. 

There exists a critical line h s in Fig. 2(a), near 
which the Meissner solution ~ 1) starts to change 
strongly. The character of this transformation depends 
on the cylinder radius R\. If the representation point 
(R\,h%) crosses the the critical line h s above the point 
Ra (R\ > Ra), the static Meissner solution becomes 
absolutely unstable relative small variations in its form, 
and transforms by a jump into new stable (static) su- 
perconducting solution with the same to. In this new 
state (which may be called as " edge-suppressed" , or " rim- 
suppressed" [25]) the order parameter ip(r) is strongly 
suppressed in some layer, which is situated near the su- 
perconductor surface, while the magnetic field penetrates 
this layer practically without screening (see Fig. 2 in [25] 
for m — and also Figs. 3-5 below). 

The jump transformation of the order parameter shape 
is acompanied by jumps in the magnetization and in 
other physical characteristics (see Fig. 3 below). The 
amplitude of these jumps diminishes with R\ diminish- 
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ing, and for R\ = Ra the jump amplitude (see the bro- 
ken vertical lines in Fig. 3) vanishes. Thus, for R\ > Ra 
the transition of the Meissner state (m = 0) to the edge- 
suppressed form (to = 0) is the first order phase transi- 
tion (by a jump). 

For R\ < Ra the Meissner solution also changes its 
form significantly, when the line h s is crossed, but this 
transformation happens smoothly, without jump (a sec- 
ond order phase transition to the edge-suppressed form). 
In the case R\ = Ra the jump amplitude is zero, so 
the first and second order phase transitions to the edge- 
suppressed state become indistinguishable. (The mean- 
ing of the points R w , i? , Ra m Figs. 2(a,d) will be 
explained later.) 

As was noted, to every point of s-region in Fig. 2(a) 
corresponds some solution of Eqs. (l)-(6) for func- 
tions ip and U, which describe the superconducting state. 
When the representation point (Rx, h^) moves, this state 
changes. Let us trace in more details, what happens, 
if the parameter R\ = R/X remains constant (i.e. the 
temperature remains constant), but the parameter fi£ — 
H/H^ changes (i.e. only the field H changes). 

In Fig. 3 are depicted as functions of H: (a) - the 
system magnetic moment (— AttMx) (8); (6) - the aver- 
age magnetic field 6 [or, the total magnetic flux inside 
the cylinder, 0i = \p^b (8)]; (c) - the difference of the 
system free energies, Ag, (10); (d) - the magnetic field 
magnitude at the cylinder axis, bo; (e) - the maximal 
value of the superconducting order parameter, ipmax', (/) 
- the order parameter value at the cylinder surface, -01- 
Fig. 3 corresponds to the case m — 0, x — 2, R\ = 5. 
[The superconducting region in Fig. 2(a) is crossed along 
the line R\ — 5.] The values of h\ = >c 2 h^ are plotted 
on the horizontal axis. By open circles the points are 
marked, which deserve a special commentary. 

As is seen from Fig. 3(a), the magnetic moment 
(— AttM\) increases linearly for small hx > (hx = 
at the point 4)- The average magnetic field inside the 
cylinder, b, is essentially reduced, in comparison to the 
normal state (see the dotted line n in Fig. 3(6)). (Analo- 
gously behaves the total magnetic flux, confined inside 
the cylinder, <j>i). The initial linear part reflects the 
Meissner effect, it corresponds to the external field ex- 
pulsion from the superconductor interior. If the field H 
is increased further, the magnetic moment increases up 
to the point 5, where the magnetization diminishes by 
jump, and takes the value 6. 

The analogous jumps are present in other quantities, 
depicted in Figs. 3(6 — /). The value of the order pa- 
rameter at the surface, -01 (Fig- 3(/)), changes espe- 
cially strongly: from ipi = 0.706 at the point 5 (where 
h x = 1.6837), to -01 = 0.087 at the point 6 (where 
h\ = 1.6838). The jump in the free energy Ag (Fig. 
3(c)) is not seen in this scale. 

When the external field fi£ is increased beyond the 
point 6, the maximal value of the order parameter, ip max , 



starts diminish quickly. This is acompanied by a fast 
drop in the "tail" of the magnetization. (This "tail" lays 
between the points 6 and 7in Fig. 3(a) and represents the 
edge-suppressed solution, eo-) When the field h% reaches 
the point 7, the order parameter vanishes completely, by 
the second order phase transition to the normal state 

(ipmax —> 0). 

Fig. 4 (to = 0, k = 2) demonstrates, what is going on, 
when the cylinder radius R is diminishing. [The super- 
conducting region in Fig. 2(a) is crossed along the line 
R\ = 3.] Now, there are no jumps on the magnetiza- 
tion curve (Fig. 4(a)), which were represented by dashed 
lines in Fig. 3(a). But, the points of inflection 2 and 4 
remain, where the magnetization curve has the maximal 
derivative. The corresponding points in Figs. 4(6, c) are 
marked by the same numerals 2 and 4 as m Fig- 4(a). 
Evidently, the solution transforms smoothly, when the 
field hx passes the critical value h s (the point 4)- This 
point corresponds approximately to the value Ra in Fig. 
2(a), where the derivative of the magnctzation curve is 
infinite. 

In Fig. 5 the case to = 0, >c = 2, Rx = 1 is presented. 
[The superconducting region in Fig. 2(a) is crossed along 
the line Rx = 1.] Here, the inflection point 4 (hx — 5.91) 
on the curve (— AttMx) practically coincides with the ter- 
mination point 5 (hx = 5.96). This situation corresponds 
to the point Rw in Fig. 2(a) (where h s and h c merge), 
and where the tail of the magnetization disappears (i.e. 
the width of the tail vanish). For Rx < R w the second 
order phase transition to the normal state occures, and 
the magnetization curve has no tail. 

Evidently, the picture, presented by Fig. 2(a) (to = 0), 
is symmetric relative the axis h^ = 0. The analogous 
picture, which corresponds to the phase sheet to = 1 
[Fig. 2(d)], is asymmetric relative the axis h% = 0. The 
edge-suppressed state exists only in the region h^ > 0, 
and they do not exist in the region h^ < 0. This may be 
interpreted in the physical terms. 

The phase sheet to = 1 describes situation, when at 
the cylinder axis there is a single vortex (its own field is 
assumed to have a positive sign), with the order param- 
eter vanishing on the axis: ipo = 0. If the external field 
of the same direction > 0) is applyed, the magnetic 
field is additionally pumped into the specimen. When 
the field reaches the critical value h s , the superconductiv- 
ity is partly destroyed, and a new superconducting state 
forms, with strongly suppressed order parameter ipi at 
the cylinder surface (see Fig. 6(/) below). This state 
preserves up to the field h^ = hi + ^ (the right branch of 
the critical line, h c > 0), until the order parameter van- 
ishes finally everywhere. Thus, in the vicinity of the line 
hi + ^ in Fig. 2(d), the mechanism of the order parameter 
suppression by the external field acts. 

If the external field has the opposite sign [h$ < 0, Fig. 
2(d)], the field of the vortex (to = 1) is gradual pumped 
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out from the vortex interior. When the external field 
reaches the critical value h% = h c ^ (the left branch of 
the critical line, h c < 0), the vortex magnetic field bo 
is completely sucked out from the superconductor (see 
Fig. 6(d)). When this happens, the vortex field at the 
axis vanishes (this means the jump transition from the 
state to = 1 to the state m = 0); the order parameter 
also vanishes by a jump (see Figs. 6(e, /)). Thus, in the 
vicinity of the line hi ^ , the mechanism of pumping out 
the field from the vortex interior acts. This explains the 
asymmetric behavior of the curves hi~ ] and h { + ] in Fig. 
2(d)- 

Note also, that the curves hi ^ and hi^ come up to 
the axis h^ = 0, having different derivatives (the fracture 
point at h^ = is marked by an open circle in Fig. 2(d); 
the corresponding value at i?A-axis is Ro). 

Mention one more peculiarity of the curves, depicted 
in Fig. 2(d). In the case of small radii (for R\ < R 
and to = 1) there exists a field interval, where the su- 
perconducting state, which was impossible at smaller h^, 
becomes possible again for larger h^. Here one has an ex- 
ample of the so-called "reentry" superconductivity. In a 
cylinder of very small radius (R\ < Ro) the magnetic 
field of the vortex can not be confined inside the su- 
perconductor and dissipates outside through the cylinder 
surface. However, the imposition of the external field pre- 
vents the dissipation of the vortex field and stabilizes the 
superconducting state. If the external field is increased 
further, the superconductivity will be finally destroyed. 
Thus, on the lowest part of the curve h c in Fig. 2(d) 
there exists a minimum at R\ = R\ m in (marked by an 
open circle). 

If the parameter x is increased, the curves h s and h c 
behave analogously to what is shown in Figs. 2(a, d). 
But, if x diminishes (see Figs. 2(6, e)), the point, where 
the curves h s and h c merge, raises up to larger R\, and 
the width of the region, where the edge-suppressed state 
exists, diminishes. At x = 1/V2 the width of this re- 
gion vanishes, and the curves h s and h c merge into one 
indivisible curve h c . For x < l/\/2 (i.e. for type-I su- 
perconductors) the edge-suppressed states do not exist. 
The behavior of the critical curves h c in the case of type- 
I superconductor is depicted in Fig. 2(c, /) for different 
x < I/V2 (to = and 1). One can see, that the super- 
conducting state with a vortex (to = I) at the cylinder 
axis, in principle, can exist in both type-I and type-II 
superconductors. 

The solution behavior in a case to = 1, x = 2 [when 
the superconducting region s in Fig. 2(d) is crossed along 
the line R\ = 5] is illustrated in Fig. 6. There are also 
jumps, as in Fig. 3, but the curves are asymmetric rela- 
tive the axis — 0. In particular, on the magnetic mo- 
ment curve (Fig. 6(a)), in a case h\ < there is no "tail", 
which exists in Fig. 3(a) for h\ < (between points 1 
and 2). This asymmetry (as was already mentioned) may 



be explained by different mechanisms, which act in de- 
struction of the superconducting state to = 1. If H > 0, 
the additional pumping of the external field into the su- 
perconductor occures, with the subsequent suppresion of 
the order parameter by the field, followed by the transi- 
tion to the normal state (or, to the state m — * m+ 1). If 
H < 0, the field is pumped out from the vortex and the 
superconductor passes into the vortex-free state to = 
(or, even into the normal state). 

The analogous picture (for m = 1, x = 2, R\ = 3) 
is presented in Fig. 7. In difference to Fig. 4 (to = 0), 
the curves in Fig. 7 (to = I) have no smooth "tails" for 
/i£ < (to the left of points 2) , but they terminate by a 
jump at the points 2, where the field at the vortex axis, 
6 , vanishes, and the order parameter suffers a jump. 

Fig. 8 corresponds to the case m = 1, x=2, R\ = 1. 
Notice, that in Fig. 8(a), at positive values h\ > 
(in the interval between the points 3 and 4) the mag- 
netic moment is positive (AttM\ > 0) , which corresponds 
to the paramagnetic susceptibility. In the interval be- 
tween the points 5 and 6 the magnetic moment is nega- 
tive (4ttM\ < 0), which corresponds to the diamagnetic 
susceptibility. As can be seen from Fig. 8(6), in the in- 
terval between points 3 and 4, the external field enhances 
the superconductivity (because ipmax increases). In the 
interval between points 4 an d 5, the external field sup- 
presses the superconductivity (because ipmax decreases). 

The analogous dependencies are presented in Fig. 9 for 
type-I superconductor (R\ = 5, x=0.5,to = 0, 1). The 
singular points are marked by an open circles. There are 
no "tails" on the magnetization curves (Fig. 9(a)) in both 
regions < and > 0. This is characteristic for type- 
I superconductors with x < l/\/2, when the destruction 
of the superconductivity by the magnetic field proceeds 
in a jump, as the first order phase transition. [This is 
valid, if R\ > Ra in Fig. 1(6). If R\ < Ra, the phase 
transition to the normal state is of second order, even in 
type-I superconductors; see Fig. 14 below.] 

There are other features, common to all type-I su- 
perconductors. Thus, if the family of curves, depicted 
in Fig. 2(e) (for to = and different x < 1/V2), is 
re-drawn on the (i?^, ft-^)-plain, one obtains the unified 
phase curve, Fig. 10(a) (to = 0). Notice also, that in the 
case to = the curves with x < l/y/2 and x > l/v2 
almost coincide, so the phase curves in Fig. 10(a) look 
like one universal function. (For > 1 this function is 
w 2.8/h$). 

Analogously, if the family of curves in Fig. 2(f ) is 
re-drawn on the (i?{, /i{)-plain, the Fig. 10(6) emerges. 
One can see from this figure, that for all x < l/y/2 there 
exists a minimal radius R^ m i n s=s 1.3. For R^ < R^ m in 
the superconducting state with a vortex (to = 1) on the 
cylinder axis is impossible. 

In Fig. 10(6) (to = 1), among others, the phase dia- 
grams are depicted for x > 1 . The corresponding curves 
(for > and R^ » 1) have the asymptote = 1 
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(or H = H C 2), common to all type-II superconductors 
with x > 1 (see Figs. 1 and 2). However, if x < 1, 
these asymptotes depend on x and began deviate from 
the value h$ = 1. [Thus, h% — > 1 for x = 1; /i£ — > 1.058 
for x = 0.9; h e -> 1.218 for x = 0.8; ft s -> 1.416 for 
x = 1/V2-] This means, that in the case to = 1 the 
value x = 1 is singular for equations (l)-(6). At this 
value the solution passes from one branch (x > 1) to an- 
other (x < 1). At the same time, the value x = l/y/2 is 
also singular, because at x = l/\/2 the second branch of 
solutions appears (or vanishes), which describe the edge- 
suppressed states [compare the behavior of the curves h 8 
in Figs. 2 (a — e)]. 

Now, we present the examples of concrete solutions of 
Eqs. (l)-(6), as functions of the space coordinate p = 
r/X (0 < r < R) for m = 1. (The solutions profiles for 
the vortex- free state (m = 0) are shown in [25], Fig. 2.) 

In Fig. 11 the solutions are shown in the case R\ = 5, 
m = 1, for several values of h\. Note the curve 3 in 
Fig. 11(a), which corresponds to h\ — and describes 
a screened Meissner-type vortex state v\ (to = 1). (In 
Fig. 6(a) to this state corresponds the point 4-) The 
curve 4 m Fig- H( a ) corresponds to the point 5 in 
Fig. 6(a), where h\ = 1.8254. Here the solution branch 
vi terminates, the temporal instability develops, and at 
h\ = 1.8255 the solution passes to new stable (static) 
branch ei (the edge-suppressed state m = 1, the point 6 
in Fig. 6(a)). The curve 5 in Fig. 11(a) corresponds to 
the value h\ — 1.8255 and describes the edge-suppressed 
state ei with a vortex (m = 1) at the cylinder axis (here 
ip(0) = and the order parameter is additionally strongly 
suppressed near the cylinder surface). 

It is expedient to trace in more details the system be- 
havior in the states with different to, when the cylinder 
radius R is fixed, and only the external field H varies. If 
the superconducting region in Fig. 1(a) is crossed along 
the line R\ = 2, the states with to = 0, 1, 2, • • • , 11 are 
possible. These states are represented in Fig. 12 by the 
following dependencies: (a) - the normalized free-energy 
difference (Ag); (6) - the magnetization (— A-kM\)\ (c) - 
the maximal value of the order parameter (ip max ) in the 
state m. 

In Fig. 12(a) the free-energy curves with the adjacent 
TO-values intersect at the points, marked by the open cir- 
cles. At these points the equilibrium transitions from the 
state to to the lower energy state m + 1 (or reversed) may 
occure. Such equilibrium transitions are acompanied by 
the reversible jumps in the magnetization (see the broken 
vertical lines in Fig. 12(6)). Thus, the equilibrium mag- 
netization curve should be reversible and have a charac- 
teristic saw-like shape. Note also, that in the case of equi- 
librium transitions the magnetization should be positive 
for all external fields (— 4nM\ > 0), this corresponds to 
the diamagnetic susceptibility. [The Meissner state with 
to = is always diamagnetic] 

However, if the metastable states with to > are also 



admissible (i.e. those parts of curves in Fig. 12(a), which 
lay to the left of the intersection points), then the mag- 
netization may turn negative within some field interval 
(— 4wM\ < for to > 1), this corresponds to the para- 
magnetic susceptibility. The magnetization curve in this 
case might be irreversible and display the hysteresis be- 
havior. 

Note also, that some of the magnetization curves in 
Fig. 12(6) (R\ = 2) have smooth "tails" (for instance, 
at to == 0) , what indicates to the presence of the edge- 
suppression effect. For larger radii (R\ > Ra ~ 4), 
instead of the inflection points, the reversible jumps in 
the magnetization may appear, which are analogous to 
those in Fig. 6(a). Probably, such qualitative details may 
be seen in the experiments with mesoscopic samples. 

As is evident from Fig. 12(c), each of the supercon- 
ducting TO-states can exist only within a limited interval 
of fields, where ip max > 0. When the field h$ is increased, 
only the last state (with to = 11) survives. In this state 
the superconducting region with ip ^ is maximally 
pushed toward the cylinder surface, but the interior is 
almost normal (see Fig. 13 below). Evidently, the state 
with to = m max is just the state of the surface supercon- 
ductivity [26,27]. In a massive cylinder (R\ 3> 1) the last 
traces of the superconductivity may persist up to the field 
h(_ = 1.69 [or, the third critical field H = H c3 = 1.69if c2 ]. 
At smaller R\ the superconductivity may persist in the 
fields, which exceed H c3 (see Fig. 1(a)). 

A peculiar field-compression effect is also present in 
Fig. 12. The first superconducting solution, which ap- 
pears in a field-cooled regime, has maximal possible to 
(the state of the surface-enhanced superconductivity). 
When the external field diminishes, this state becomes 
impossible, but new to — 1-state appears, in which the su- 
perconducting region is shifted toward the center, and si- 
multaneously the edge-suppressed region begins to form. 
The state to = 4, for instance, has both the region of 
edge-suppressed superconductivity (where the external 
field is practically not screened), and the region at some 
distance from the center, where the field of the giant vor- 
tex (to = 4) is compressed and the superconductivity 
is also almost crushed (because at the center the order 
parameter is always small, ip(p) ~ p m <C 1). (A flux- 
compression effect was discussed also in [6]). 

As was mentioned, Fig. 12(6) indicates to the possible 
paramagnetic susceptibility of the superconducting cylin- 
der in the state with to > 0. This topic also deserves a 
special comment. 

The paramagnetic Meissner effect (or, Wohlleben ef- 
fect) was initially observed in [40] and extensively 
discussed in the literature (see, for example, Refs. 
[6,8,11,15,41-43] and references therein). A number of 
controversial explanations of this effect were proposed, 
but there is as yet no clear understanding of its nature. 
This effect is observed, in particular, in a finite-dimension 
samples, in the field-cooled regime, it shows the signs of 
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metastability, of the reentrant behavior and hysteresis 
[43]. It is interesting, that all these features follow natu- 
rally from Eqs. (l)-(6) of the Ginzburg-Landau theory. 

Consider, for instance, the state to = 4 (which is dis- 
tinguished, for clarity, by thick curves in Fig. 12). This 
state can not exist in small external fields fi£, because 
the vortex magnetic field (m = 4, R\ = 2) is too strong 
for a mesoscopic sample, and can not be confined inside. 
However, the imposition of a finite external field prevents 
the field dissipation, the situation stabilizes and the state 
to = 4 becomes possible again in a larger field (the reen- 
trance effect). Evidently, it is more easy to occupy the 
state to = 4, going down from the stronger fields (in a 
field-cooled regime) . The negative part of the magnetiza- 
tion curve to = 4 (the paramagnetic branch) corresponds 
to a metastable state of higher free-energy (see Fig. 12. 
(a)), than the state m = 3. [Here the state m = 4 is de- 
fined as metastable, on the ground, that its free energy 
is less than the free energy of the state to = 3. Such 
metastability means, that the irreversible transition to 
the lower laying state to — 3 is possible. However, the 
state m — 4 continues to be stable relative small pertur- 
bations of its form.] 

Note, that in the positive part of the magnetiza- 
tion curve the superconducting currents flow, mainly, to 
screen out the external field. [This leads to the diamag- 
netic response in Fig. 12(6).] In the negative part of the 
magnetization curve the role of the vortex field is more 
important, and the superconducting screening currents 
flow, mainly, in the opposite direction. [This means the 
sign reversion of the magnetic moment, because it is pro- 
portional to the current, M ~ f[j s r]dv. The presence of 
positive and negative contributions to the total magnetic 
moment is evident also from Fig. 11(/), where the cur- 
rent density may be positive or negative). At the point 
M\ = the currents, which screen the internal and exter- 
nal fields, counterbalance each other. Thus, according to 
the Ginzburg-Landau theory, the paramagnetism present 
in Fig. 12, has purely electrodynamic explanation. 

Now we shall demonstrate, how the order parameter 
vanishes in approaching the critical curves h!f^ > in 
Figs. 1 and 2 [the curves h c are supplied with an addi- 
tional index [m)\. The order parameter is small in the 
vicinity of h[ m \ In Fig. 13 the degenerated solutions 
(i.e. those with ip max <C 1) are shown for R\ = 2, x = 2 
and different to. [The curves in Fig. 13 are normalized 
to unity, with true values ipmax « lj In the origin of 
coordinates the order parameter behaves as ip(p) <~ p m , 
in agreement with the theory (see, for instance, [29,35]). 
At p > 1 the order parameter is also small (here the 
edge-suppressed state e m is formed) . For all degenerated 
solutions, depicted in Fig. 13 (with ip max "C 1), the field 
B(r) w H. In this case Eqs. (1),(2) can be linearized 
and reduced to one linear equation for the function tjj 
[26,27]. The solutions of this equation may be found 



analitically, in terms of the Kummer functions (see, for 
instance, [6,8]). [Notice, however, that the normalization 
factor in the linear equation for tp remains arbitrary, and 
to find the true value ip max it is necessary to solve the 
full system of nonlinear equations (l)-(6).] The critical 
fields hi m ^ > in Fig. 1, which we found by solving the 
full system of equations (1)— (6), may be obtained more 
easily on the base of a linear equation (when tp m ax <C 1), 
as was done in [6,8]. [However, when approaching the 
phase boundaries h[ m ^ for negative < 0, or to find 

the critical fields h{ m \ when the edge-suppressed state 
forms (see Figs. 1 and 2), but the order parameter ip max 
remains finite, it is necessary to use the full system of 
equations (l)-(6).] 

It is evident from Fig. 13, that even the degener- 
ated solutions (with ip max <C 1) display the presence 
of the edge-suppresssion effect. If m = m max , the or- 
der parameter is maximal at the cylinder surface. For 
m < m max the order parameter begins being suppressed 
at the boundary. When the field H diminishes, only the 
states with relatively small to survive [with the inflexion 
points on the magnitization curves, see Figs. 12(6, c)], 
and finally the edge-suppressed solutions would form 
[analogous to those, shown in Fig. 11, with the char- 
acteristic jumps, as in Figs. 3 and 6]. 

Finally, it is of interest to note, that the transitions 
between states with different to are also possible in type-I 
superconductors (see also [35,11]). This is demonstrated 
by Fig. 14. In Figs. 14(a, 6) the profiles of the self- 
consistent solutions ip{r) and 6(r) are shown for type-I 
superconductor (x = 0.5, R\ = 4) in the fields: — 
1.4892 (to = 0); h 6 = 2.0344 (to = 1); h 6 = 2.1084 
(to = 2). These ft-j-values are in the immediate vicinity 
of the corresponding critical curves hi m) in Fig. 1(6): if 
/i£ is increased by 1 • 10~ 4 [i.e., if the critical curves in Fig. 
1(6) are crossed along the line R\ = 4], the transitions 
to the normal state occures (ip(r) = 0, B(r) = H). 

As is clear from Fig. 14(a), in type-I superconductors 
the transition to the normal state is of the first order (by 
a jump from a finite value i>max)- However, the jump 
amplitude (i.e. the value ipmax a ^ t ne transition point) 
depends on parameters m, x and R\ (or R^ — R\/x). 
The jump amplitude, ^^ ax {R\), is depicted in Fig. 14(c) 
for k = 0.1 and k = 1/V2 in the case to = 0. As follows 
from Fig. 14(c), the value ip^ax diminishes, when R\ di- 
minishes, with the jump amplitude vanishing at the point 
Ra- Thus, the points Ra on the critical curves in Figs. 
1(6), 10(a) and 14(c) divide two regions of the phase tran- 
sitions in the case of type-I superconductors: for larger 
R\ > Ra the first order phase transition to the normal 
state occures; for smaller R\ < Ra one has the second 
order phase transition. At R\ = Ra both transitions 
became indistinguishable, because the jump amplitude 
ipm ax vanish. [Remind, that in type-II superconductors, 
with x > l/\/2, the phase transition to the normal state 
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is always second order, regardless of the specimen size.] 
The critical size Ra, which separates the regions of the 
first and second order phase transitions in the case of 
type-I superconductors, may also be found from different 
considerations [1,32,36,38]. Fig. 14(c) shows, that this 
size depends slightly on x. 

IV. COMMENTS AND CONCLUSIONS 

In conclusion, the main results of the present investi- 
gation are formulated. 

The self-consistent solutions of a nonlinear system of 
Ginzburg-Landau equations, which correspond to the 
fixed values of the magnetic quantum number m and have 
a cylinder symmetry (i.e. depend only on the radial coor- 
dinat r, but not on the azymuthal angle ip), are studied. 

It is found, that for type-II superconducting cylinder 
of finite radius (R > Ra), placed in sufficiently weak 
magnetic field H, the order parameter ip(r) ~ const near 
the cylinder surface, and behaves as ip(r) ~ r m near the 
center (the giant- vortex state). If the field is increased 
(H > the giant- vortex state changes abruptly 

and acquires new form. (A first order phase transi- 
tion to the edge-suppressed state occures.) In this new 
state the order parameter is strongly suppressed near the 
cylinder edge, with the magnetic field penetrating prac- 
tically without screening into some layer in the vicinity 
of the surface. At the center the behavior ij)(r) <~ r m 
remains. This transition to the edge-suppressed form is 
acompanied by jumps in the magnetization and in the 
total magnetic flux, confined in the cylinder. The prop- 
erties of such edge-suppressed solutions (which depend 
on the parameters m, H, x, T) are studied. It is shown, 
that the edge-suppressed states exist only in type-II su- 
perconductors, but the giant-vortex states are possible 
both in type-I and type-II superconductors. The phase 
boundaries are found, inside which the giant m- vortex so- 
lutions can exist. In particular, the characteristic radius 
of the cylinder is found, Ra, which separates the regions 
of the first and second order phase transitions between 
the superconducting and normal states. 

The behavior of the magnetization (M) and of the to- 
tal magnetic flux (<I>i) as functions of the external field 
(H) is found, with account of the jumps, which acom- 
pany the formation of the edge-suppressed states. It is 
demonstrated, that some maximal value m = m max ex- 
ist, for which the edge-suppressed state degenerates into 
the usual state of the surface superconductivity (with 
nucleation field H sa H c3 [26,27]). The Gibbs free en- 
ergies of different axially symmetrical superconducting 
states are compared, and the problem of paramagnetic 
effect, observed in mesoscopic samples, is shortly dis- 
cussed. It is found, that the paramagnetic effect is possi- 
ble in metastable vortex states and it may be attributed 



to the disbalance of the superconducting currents, which 
screen the external field H and the field of the vortex 
from entering the sample interior. 

As was pointed out earlier [25], according to the 
Ginzburg Landau theory, two competing mechanisms of 
the external field penetration into the superconducting 
cylinder exist - in a form of vortices [2-15,30], and in a 
form of the edge-suppressed layer. The question, which 
of these mechanisms is more favorable energetically, may 
be answered only after the full comparison is made of the 
Gibbs free energies for the arbitrary multy- vortex config- 
urations in a finite radius cylinder. This difficult the- 
oretical problem (as well as the generalisation to other 
geometries) is left outside the scope of the present inves- 
tigation. Some of the edge-suppressed states, described 
above, are metastable, but they may be observable in a 
specially arranged experiments. We beleve it is expedient 
to draw attention to this possibility. 
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Figures captions 

Fig. 1. (a) - The phase diagram of type-II super- 
conductor (x = 2). The numerals at the curves - the 
values of the magnetic quantum number to. The value 

= 1 corresponds to the critical field H c2 = 0o/(27r£ 2 ). 
The value = 1.69 corresponds to the maximal field 
H c3 = 1.69H c2 , in which the surface superconductivity is 
still possible in macroscopic specimens, (b) - The phase 
diagram of type-I superconductor (x = 0.5, m = 0, 1,2). 
The arrows and the letter i?A at the curves denote the 
points, where the first and second order phase transi- 
tions to the normal state became idistinguishable (see 
the text). 



Fig. 2. (a) - The case to = 0, x = 2. The phase 
curve h c divides superconducting (s) and normal (n) 
states of the cylinder. When the representation point 
(Rx, h^) crosses the line h s (the dotted line), a new edge- 
suppressed solution forms, with suppressed order param- 
eter -01 (see Fig. 3(/)). No edge-suppressed states exist 
for R\ < R w . (d) - Analogously, for the case to = 1, 
x = 2. When R\ < R Q , the field h$ > stimulates the 
superconducting state, so the reentry superconductivity 
is possible, (b) and (e) - Analogously for different x > 1. 
(c) and (/) - The phase curves of type-I superconductor 
(for various x < 1/V2). The edge-suppressed states in 
type-I superconductors are abscent. 

Fig. 3. The case to = 0, x = 2, R x = 5. (a) - 
The cylinder magnetic moment, M\ = M/H\. (b) -The 
average magnetic field in the cylinder, b (or, the total 
magnetic flux in the system, (f>\ = \p^b (8)). (c) - The 
normalized difference of free energies, Ag, (10). (d) - 
The magnetic field at the cylinder axis, bo. (e) - The 
maximal value of the order parameter, tp ma x- (/) - The 
order parameter at the cylinder surface, tpi. The letter 

denotes the Meissner state (vortex- free, to = 0). The 
letter eo denotes the edge-suppressed state (vortex-free, 
to = 0). The transition from one branch of the solution 
to another takes place between points 2,3, or 5,6. (The 
solution profiles for the case R\ = 5, to = 0, x = 2 are 
depicted in [17], Fig.2.) 

Fig. 4. Analogous to Fig. 3, but for the case to = 0, 
x = 2, R x = 3. 

Fig. 5. Analogous to Fig. 3, but for the case to = 0, 
K= 2, R x = 1. 

Fig. 6. The same, as in Fig. 3, but for to = 1, x = 2, 
Rx = 5. The letter vi denotes the Meissner- type vortex 
state (with a vortex to = 1 at the cylinder axis) . The let- 
ter ei denotes the edge-suppressed state with the vortex 
(to = 1) at the axis. (The solutions profiles in the case 
to = 1, x =2, Rx = 5, are shown below in Fig. 11.) 

Fig. 7. Analogous to Fig. 3, but for the case to = 1, 
x= 2, Rx = 3. 

Fig. 8. Analogous to Fig. 3, but for the case to = 0, 
x= 2, R x = 1. 

Fig. 9. The same, as in Fig. 3, but for type-I super- 
conductor (x = 0.5, Rx = 5, to = and to = 1). There 
are no edge-suppressed states and no "tails" in the mag- 
netization for type-I superconductors. 

Fig. 10. (o) - The phase diagram for to = and dif- 
ferent x. The arrow and letter i?A devide the regions of 
the first and second order phase transitions to the normal 
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state in the case of type-I superconductors. For m = 
the curves with different x [see numbers in Fig. 10(a)] 
practically coincide. (6) - The phase diagram for m = 1 
and different x. For R^ < R^ m in the superconducting 
state with the vortex (m = 1) at the axis is impossible. 

Fig. 11. The coordinate dependence of solutions for 
m = 1, R\ = 5, x = 2 and different fields (sec the values 
in Fig. 11(a). (a) - The order parameter, ip(p); p = r/X. 
(6) - The derivative, dip /dp. (c) - The second deriva- 
tive, d 2 ip/dp 2 . (d) - The field potential, U(p). (e) - The 
normalized magnetic field, b(p). (/) - The normalized 
current, j(p). 

Fig. 12. (a) - The free energy, Ag; (6) - the magnetic 
moment, (— 4nM\) and (c) - the maximal value of the 
order parameter, ipi, for the cylinder with R\ = 2, x = 2 
and various m, as functions of /ij. 

Fig. 13. The order parameter profiles in a case of de- 
generated states {ipmax <C 1) for R\ = 2, x = 2 and 
various m. (All the curves represent the self-consistent 
solutions of Eqs. (l)-(6), but they are normalized to 
unity.) With m increasing, the value i\) max approaches 
the boundary r = R. The last of the degenerated states 
(m = 11), which is maximally pushed toward the bound- 
ary, corresponds to the state of the surface superconduc- 
tivity [26,27]. 

Fig. 14. (a) - The order parameter profiles in a case of 
type-I superconductor (x = 0.5, R/X = 4). The external 
fields (see the numerals at the curves) are in the imme- 
diate vicinity of the critical boundaries hi™ 1 ^ in Fig. 1 (6) . 
If the field is increased by 1 • 10~ 4 , the jump transi- 
tion to the normal state occures. (b) - The corresponding 
magnetic field profiles, (c) - The jump amplitude, tpmax> 
as function of R\ for x = 0.1 and x = 1/V2. [The jump 
amplitude V'maz i s tnc va hic of ip max m a field, which 
just preceeds the first order phase transition to the nor- 
mal state. The value Ra, where the jump amplitude 
vanish (see Figs. 1(6), 10(a), 14(c)) depends on x.] 
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